Mon. Not. R. Astron. Soc. 296, 206-216 (1998) Printed 1 February 2008 (MN MfeK style file vl.4) 



Evolution of DA white dwarfs in the context of a new 
theory of convection 

L. G. Althaus* and O. G. Benvenutot 

Facultad de Ciencias Astronomicas y Geofisicas, Paseo del Bosque S/N, (1900) La Plata, Argentina 
Accepted 1997 November 18. Received 1997 November 18; in original form 1997 July 8 

ABSTRACT 

In this study we compute the structure and evolution of carbon - oxygen DA (hydrogen 
- rich envelopes) white dwarf models by means of a detailed and updated evolutionary 
code. We consider models with masses from 0.5 to 1.0 Mq and we vary the hydrogen 
layer mass in the interval of 10~ 13 < Mh/M < 10~ 4 . In particular, we treat the energy 
transport by convection within the formalism of the full - spectrum turbulence theory, 
as given by the Canuto, Goldman & Mazzitelli (CGM) model. We explore the effect of 
various hydrogen layer masses on both the surface gravity and the hydrogen burning. 
Convective mixing at low luminosities is also considered. 

One of our main interests in this work has been to study the evolution of ZZ 
Ceti models, with the aim of comparing the CGM and mixing - length theory (MLT) 
predictions. In this connection, we find that the temperature profile given by the CGM 
model is markedly different from that of the ML1 and ML2 versions of the MLT. In 
addition, the evolving outer convection zone behaves differently in both theories. 

We have also computed approximate effective temperatures for the theoretical 
blue edge of DA instability strip by using thermal - timescale arguments for our 
evolving DA models. In this context, we found that the CGM theory leads to blue 
edges that are cooler than the observed ones. However, because the determination 
of atmospheric parameters of ZZ Ceti stars is dependent on the assumed convection 
description in model atmosphere calculations, observed blue edges based on model 
atmospheres computed considering the CGM theory arc required in order to perform 
a self - consistent comparison of our results with observations. Finally, detailed non - 
adiabatic pulsational computations of ZZ Ceti models considering the CGM convection 
would be necessary to place the results found in this paper on a firmer basis. 
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1 INTRODUCTION 

White dwarf stars with hydrogen lines (spectral type DA) 
constitute the most numerous group of the observed white 
dwarfs. Since the first DA white dwarf was reported to 
exhibit multiperiodic luminosity variations (see McGraw 
1979), the interest in these objects has greatly increased. 
Rapid progress in this field has been possible thanks to the 
development of powerful theoretical tools as well as to a in- 
creasing degree of sophistication in observational techniques. 
Over the years, various studies have presented strong evi- 
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dence that pulsating DA white dwarfs (DAV) or ZZ Ceti 
stars represent an evolutionary stage in the cooling history 
of the majority of, if not all, DA white dwarfs (Weidemann 
& Koester 1984 and references therein). 

Studies carried out notably by Dolez & Vauclair (1981) 
and Winget et al. (1982) on the basis of detailed non - radial, 
non - adiabatic pulsation calculations, demonstrated that 
the K-mechanism operating in the hydrogen partial ioniza- 
tion zone is responsible for the g(gravity)-mode instabilities 
in ZZ Ceti stars. From then on, pulsating white dwarfs have 
captured the interest of numerous investigators, who have 
employed the powerful technique of white dwarf seismology 
to derive fundamental parameters of these stars, such as 
the stellar mass, chemical composition and stratification of 
the outer layers (see, for example, Winget et al. 1994 and 
Bradley 1996 in the context of DB and DA, respectively). In 
particular, pulsation calculations have shown that the theo- 
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retical determination of the location of the blue (hot) edge 
of the instability strip (which marks the onset of pulsations) 
in the Hertzsprung - Russell diagram is mostly sensitive to 
the treatment of convection in envelope models. Accordingly, 
a comparison with the location of the hottest observed ZZ 
Ceti stars allows us to obtain valuable information about 
the efficiency of convection in such stars (Winget et al. 1982; 
Tassoul, Fontaine & Winget 1990; Bradley & Winget 1994; 
Bradley 1996 and references cited therein). 

Unfortunately, the determination of atmospheric pa- 
rameters of ZZ Ceti stars is strongly dependent on the as- 
sumed convection description in model atmosphere calcu- 
lations, against which the observed spectrum is compared. 
The theoretical analysis carried out by Bergeron, Wesemael 
& Fontaine (1992a) is particularly noteworthy in this regard. 
In fact, these authors demonstrated that the predicted emer- 
gent fluxes, color indices and equivalent widths of DA white 
dwarfs, notably of the ZZ Ceti stars, are very sensitive to the 
details of the parametrization of the mixing - length theory 
(MLT) of convection. In particular, Bergeron et al. (1992a) 
concluded that the uncertainties in the observational defini- 
tion of the blue edge of the ZZ Ceti instability strip brought 
about by different convective efficiencies may be appreciable, 
irrespective of the observational technique that is used. 

Needless to say, the employment of MLT to deal with 
convection represents a serious shortcoming in theoretical 
studies of white dwarfs. In particular, MLT approximates 
the whole spectrum of eddies necessary to describe the stel- 
lar interiors by a single large eddy. This representation intro- 
duces in the description of the model certain free parameters 
(Bohm - Vitense 1958) which, in most stellar studies, are re- 
duced to a single one: the distance I travelled by the eddy 
before thermalizing with the surrounding medium. I is writ- 
ten as I = aH p where H p is the pressure scaleheight and 
a is the free parameter. Obviously, the existence of such an 
adjustable parameter prevents stellar evolution calculations 
from being completely predictive. In white dwarf studies, 
several parametrizations of the MLT are extensively em- 
ployed. Specifically, the ML1, ML2 and ML3 versions, which 
differ in the choice of the coefficients appearing in the expres- 
sions of the convective flux, average velocity and convective 
efficiency (see Tassoul et al. 1990 for details) are associated, 
respectively, with increasing convective efficiency (for the 
ML1 and ML2 versions a = 1, while ML3 is the same as 
ML2 but with a = 2). 

Over the last two decades, an unparalleled effort has 
been devoted to the observational determination of the lo- 
cation of the blue edge of the instability strip of the ZZ 
Ceti stars (see Wesemael et al. 1991 for a review). Based on 
different observational techniques, the majority of the stud- 
ies cited by Wesemael et al. point towards a temperature 
of the blue edge, defined by the ZZ Ceti star G117-B15A, 
of approximately 13000 K. New observational data and im- 
proved model atmospheres tend to indicate a significantly 
lower temperature for the blue edge. In fact, from ultravio- 
let observations performed with the Hubble Space Telescope, 
Koester, Allard & Vauclair (1994) have shown that model 
atmospheres calculated with the ML1 parametrization do 
not fit satisfactorily to the ultraviolet spectrum of G117- 
B15A. Instead, a more efficient convection, as given by the 
ML2 version, is required to explain the observed spectrum. 
On the basis of this result, Koester et al. derived for G117- 



B15A an effective temperature of 12250 ± 125 K. In view of 
the fact that the ultraviolet spectra of other candidates for 
the blue edge of the ZZ Ceti instability strip are similar to 
that of G117-B15A, Koester et al. concluded that the blue 
edge temperature is likely to be around 12250 K, a value 
which is considerably lower than found in previous studies. 

In a still more recent analysis, Bergeron et al. (1995) in- 
ferred an even lower effective temperature for G117-B15A. 
More precisely, on the basis of new optical and ultravio- 
let analyses of ZZ Ceti stars, Bergeron et al. succeeded in 
constraining the convective efficiency used in model atmo- 
sphere calculations of these stars. Indeed, they demonstrated 
that a parametrization less efficient than ML2, particularly 
the ML2/a=0.6 version (intermediate in efficiency between 
ML1 and ML2) yields ultraviolet temperatures that are com- 
pletely consistent with the optical determinations. In this 
respect, they derived for G117-B15A an effective tempera- 
ture of 11620 K. In the sample analysed by Bergeron et al., 
the blue edge of the instability strip is defined by the star 
G226-29, for which they obtained a temperature of 12460 
K and a surface gravity of loggr = 8.29. Finally, Fontaine 
et al. (1996) presented a reanalysis of the observed pulsa- 
tion modes of G117-B15A that reinforces the conclusions 
of Bergeron et al. (1995) on the atmospheric parameters of 
this variable white dwarf. Fontaine et al. concluded that a 
model atmosphere of G117-B15A characterized by an effec- 
tive temperature of « 11500 K and a ML2/a=0.6 convection 
is consistent with optical and ultraviolet spectroscopic ob- 
servations as well as with the observed pulsation amplitudes 
in different wavelength bands. 

Another conclusion drawn by Bergeron et al. (1995) (see 
also Wesemael et al. 1991) worthy of comment is related 
to the fact that convective efficiency in the upper atmo- 
sphere appears to be quite different from that in the deeper 
envelope. In fact, according to non - adiabatic pulsation 
calculations (Bradley & Winget 1994 and references cited 
therein), a more efficient convection than that provided by 
ML2/a=0.6 is needed in the deeper envelope to explain the 
observed DA blue edge. This result, which reflects the in- 
ability of MLT to describe convection throughout the entire 
outer convection zone adequately, has indeed been borne out 
by recent detailed hydrodynamical simulations of convection 
in ZZ Ceti stars (see later in this section). Other objections 
have been raised against MLT. For instance, the value of 
q, which is usually derived from solar radius adjustments, 
is larger than unity, in contrast to the basic postulates of 
MLT. In addition to this, the procedure of applying the so- 
lar a parameter to other stars is not entirely satisfactory, at 
least in the context of red giant stars for which a wide range 
of a values is needed (Stothers & Chin 1995, 1997). Finally, 
MLT does not fit laboratory convection data (Canuto 1996). 

Fortunately, a considerable effort has been devoted in 
recent years to improving some of the basic postulates of 
MLT. In particular, Canuto & Mazzitelli (1991, 1992, here- 
after CM) and more recently, Canuto, Goldman & Mazz- 
itelli (1996, hereafter CGM) developed two convection mod- 
els based on the full - spectrum turbulence (FST) theory, 
which represent a substantial progress in modelling stellar 
convection (we refer the reader to Canuto & Christensen - 
Dalsgaard 1997 for a recent review on local and non - local 
convection models). Both of the treatments resolve essen- 
tially two main shortcomings of the MLT 
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(i) The MLT one - eddy approach is replaced by a full 
spectrum of turbulent eddies derived on the basis of amply 
tested, modern theories of turbulence. As discussed by CM, 
the one - eddy approximation is only justified for viscous flu- 
ids but becomes completely inadequate for almost - inviscid 
stellar interiors. 

(ii) There are no adjustable parameters. 

These features render the CM and CGM models sub- 
stantially more solid and complete than MLT. The CM 
model in particular has successfully passed a wide variety of 
laboratory and astrophysical tests, performing in all cases 
much better than MLT. With regard to the astrophysical 
case, the CM model has been applied to the study of dif- 
ferent kind of stars and evolutionary phases. Applications 
include the solar model (see CM), the effective tempera- 
ture of which is fitted within w 0.5 per cent accuracy; red 
giants (D'Antona, Mazzitelli & Gratton 1992 and Stothers 
& Chin 1995, 1997), for which MLT may require various 
values of a for different stars, helioseismology (Paterno et 
al. 1993 and Monteiro, Christensen-Dalsgaard & Thompson 
1996 amongst others), low - mass stars (D'Antona & Mazz- 
itelli 1994) and stellar atmospheres (Kupka 1996). In the 
white dwarf context, the CM model has been shown to be 
a very valuable tool as well. In fact, Althaus & Benvenuto 
(1996) (see also Mazzitelli & D'Antona 1991) studied the 
evolution of carbon - oxygen white dwarfs with helium - rich 
envelopes (spectral type DB) and found that the CM model 
predicts theoretical blue edges for the DB instability strip 
in good agreement with the observations of Thejll, Vennes 
& Shipman (1991). Mazzitelli (1995) has also applied the 
CM model to the study of convection in DA white dwarfs. 
More recently, Althaus & Benvenuto (1997a) used the CM 
model to study the evolution of helium white dwarfs of low 
and intermediate masses. 

CGM have recently improved the CM model by develop- 
ing a self-consistent model for stellar turbulent convection, 
which, like the CM model, includes the full spectrum of ed- 
dies but computes the rate of energy input self - consistently. 
More precisely, CGM resolve analytically a full turbulence 
model in the local limit, taking into account the fact that 
the energy input from the source (buoyancy) into the turbu- 
lence depends on both the source and the turbulence itself. 
At low and intermediate convective efficiencies, the CGM 
gives rise to higher convective fluxes than those given by 
the CM model. CGM find that the main sequence evolu- 
tion of a solar model does not differ appreciably from the 
CM predictions but demands a smaller overshooting to fit 
the solar radius, in agreement with new observational data. 
D'Antona, Caloi & Mazzitelli (1997) have also applied the 
CGM model to study the problem of globular cluster ages. 
In the white dwarf domain, the CGM model also fares bet- 
ter than the CM model (Althaus & Benvenuto 1997b and 
Benvenuto & Althaus 1997). Finally, Benvenuto & Althaus 
(1997b) have applied the CGM model to study the evolution 
of low - and intermediate - mass helium white dwarfs with 
hydrogen envelopes. 

In order to assess to what extent the MLT is able to 
provide a reliable description of convection in the outer lay- 
ers of DA white dwarfs, Ludwig, Jordan & Steffen (1994) 
have performed detailed two - dimensional hydrodynamical 
simulations of convection at the surface of a ZZ Ceti star. 
They demonstrated that the temperature profile in the outer 



layers of a T ff = 12600 K, log g = 8.0 DA model cannot be 
reproduced with model atmospheres calculated with MLT 
and with a single value of a. Instead, the value of a is to 
be varied from 1.5 in the upper atmosphere to 4 in much 
deeper layers in order for the hydrodynamical temperature 
profile to be reasonably well represented. Another result ob- 
tained by Ludwig et al. is that their hydrodynamical sim- 
ulations lead to convective zones that are too shallow for 
the k, 7 mechanism to give rise unstable g modes, imply- 
ing a blue edge for the ZZ Ceti instability strip significantly 
cooler than 12600 K. This important finding has recently 
been borne out by Gautschy, Ludwig & Freytag (1996) on 
the basis of non - radial, non - adiabatic pulsation calcula- 
tions. Gautschy et al. included in their analysis Ludwig et 
al.'s hydrodynamical simulations to derive the structure of 
the convective layers of their white dwarf models. This fea- 
ture represents substantial progress compared with earlier 
investigations, which employ MLT to deal with convection. 
Gautschy et al. are unable to find a theoretical blue edge 
for the ZZ Ceti instability strip consistent with the observa- 
tions of Bergeron et al. (1995) and conclude that, according 
to their calculations, it lies between 11400 K and 11800 K 
(for a model with log g = 8). 

The present study is devoted to presenting new calcu- 
lations of the evolution of carbon - oxygen DA white dwarfs 
considering the model for stellar turbulent convection of 
CGM as well as the various parametrizations of MLT com- 
monly used in the study of these stars. Attention is focused 
mainly on those evolutionary stages corresponding to the 
ZZ Ceti effective temperature range. The calculations arc 
carried out with a detailed white dwarf evolutionary code in 
which we include updated equations of state, opacities and 
neutrino emission rates. Furthermore, we take into account 
crystallization, convective mixing and hydrogen burning. To 
explore the sensitivity of the results to various input model 
parameters, we vary both the model mass from 0.50 - 1.0 
Mq , which covers most of the observed mass distribution of 
DA white dwarfs (Bergeron, Saffcr & Licbcrt 1992b, Marsh 
ct al. 1997), and the hydrogen layer mass Mh in the inter- 
val 10~ 13 < Mh/M < 10~ 4 where M is the model mass. 
Finally the hydrogen - helium transition zone is assumed 
to be almost discontinues. Details of our evolutionary code 
and its main physical constituents are presented in Section 
2. Section 3 is devoted to analysing our results and, finally, 
in Section 4 we summarize our findings. 



2 COMPUTATIONAL DETAILS AND INPUT 
PHYSICS 

The white dwarf evolutionary code we employed in this 
study was developed by us independently of other re- 
searchers and it has been used to study different problems 
connected to white dwarf evolution as well as the cooling 
of the so - called strange dwarf stars (Benvenuto & Althaus 
1996). Details of the code and its constitutive physics can 
be found in Althaus & Benvenuto (1997a) and Benvenuto & 
Althaus (1997, 1998). In what follows we restrict ourselves 
to a few brief comments. 

The equation of state for the low - density regime is 
that of Saumon, Chabrier & Van Horn (1995) for hydrogen 
and helium plasmas. The treatment for the high - density, 
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completely ionized regime appropriate for the white dwarf 
interior includes ionic and photon contributions, coulomb 
interactions, partially degenerate electrons and electron ex- 
change and Thomas - Fermi contributions at finite tempera- 
ture. Radiative opacitites for the high - temperature regime 
(T > 8000 K) are those of OPAL (Iglesias & Rogers 1993) 
with metallicity Z — 0.001, whilst for lower temperatures 
we use the Alexander & Ferguson (1994) molecular opaci- 
ties. Molecular effects become relevant at temperatures as 
high as 5000 K and below 2500 K they are dominat (Alexan- 
der & Ferguson 1994). We shall see that convective mixing 
at low effective temperatures may in some cases change the 
outer layer composition from almost pure hydrogen to al- 
most pure helium. This being the case, we had to rely, in 
the low - temperature regime, on the older tabulation of Cox 
& Stewart (1970) for helium composition. Conductive opac- 
ities for the low - density regime are taken from Hubbard 
& Lampe (1969) as fitted by Fontaine & Van Horn (1976). 
Conductive opacities for the liquid and crystalline phases, 
and the various mechanisms of neutrino emission relevant 
to white dwarf interiors are taken from the works of Itoh 
and collaborators (see Althaus & Benvenuto 1997a and ref- 
erences cited therein). It is worth mentioning that neutrino 
cooling is dominant during the hot phases of white dwarf 
evolution but they become completely negligible at the ZZ 
Ceti effective temperature range. 

The procedure we followed to generate the initial mod- 
els of different masses is described in Benvenuto & Althaus 
(1995). Such initial models were constructed starting from 
a 0.55Mq carbon - oxygen white dwarf model kindly pro- 
vided to us by Professor Francesca D'Antona. We assumed 
the same core chemical composition (consisting of a mixture 
of carbon and oxygen, see Fig. 1 of Benvenuto & Althaus 
1995) for all of them, notwithstanding the changes in the in- 
ternal chemical composition that are expected to occur for 
models with different masses because of differences in the 
pre - white dwarf evolution of progenitors objects. The car- 
bon - oxygen core of our models is surrounded by an almost 
pure helium envelope, the mass of which (Mh 8 ) was varied 
in the range 10~ 6 < M Ho /M < 10~ 2 . In DA white dwarfs 
there is an almost pure hydrogen envelope on the top of 
the helium layers. Unfortunately, the mass of this hydrogen 
envelope is only weakly constrained by pre - white dwarf 
evolutionary calculations (DAntona & Mazzitelli 1991) and 
by non - adiabatic pulsation studies (Fontaine et al. 1994). 
In recent years, however, strong evidence has been accumu- 
lated in favour of the idea that some ZZ Ceti stars appear 
to have thick hydrogen layers (Fontaine et al. 1994). In the 
present study, we decided to treat the mass of the hydro- 
gen envelope as basically a free parameter within the range 
10~ 13 < M H /M < 10~ 4 . We refer the reader to Benvenuto 
& Althaus (1998) for details about the procedure we fol- 
low to add a hydrogen envelope at the top of our models. 
We want to mention that we assumed the hydrogen/helium 
transition zone to be almost discontinuous. 

We also included in our evolutionary calculations the 
release of latent heat during crystallization (see Benvenuto 
& Althaus 1997) and the effect of the convective mixing. 
We shall see in particular that the onset of crystallization 
in our more massive models occurs at effective temperatures 
close to the observed blue edge. Finally, we also considered 
the effects of hydrogen burning by including in our numer- 



ical code the complete network of thermonuclear reaction 
rates for hydrogen burning corresponding to the proton - 
proton chain and the CNO bi - cycle. Nuclear reaction rates 
are taken from Caughlan & Fowler (1988) and j3 - decay 
rates from Wagoner (1969), taking into account the correc- 
tions for their Q - values resulting from the effect of neu- 
trino losses. Electron screening is from Wallace, Woosley & 
Weaver (1982). We consider the following chemical species: 
X H, 2 H, 3 He, 4 He, 7 Li, 7 Be, 8 B, 12 C, 13 C, 13 N, 14 N, 15 N, 
15 0, ie O, 17 and 17 F. 

The most relevant feature of the present study is the em- 
ployment of the FST theory to deal with energy transport 
by convection. This theory represents a substantial improve- 
ment over MLT, particularly in the treatment of low - viscos- 
ity fluids like the ones present in stellar interiors, for which 
the MLT one - eddy assumption is completely inadequate. 
In particular, a self - consistent model without adjustable 
parameters based on the FST approach has been recently 
formulated by CGM. These authors fitted their theoretical 
results for the convective flux F c with the expression 



F c = XT//- 1 (V-V ad )$, 



(1) 



where K — AacT 3 /3Kp is the radiative conductivity, H p is 
the pressure scaleheight, and V and V a d are, respectively, 
the true and adiabatic temperature gradients. $ is given by 



$ = F!(S)F 2 (S), 
where 

Kn 

Fi(S) = (^) 3 aS k [(l + bS) 1 -l] q 
l.o 

and 

c.0.72 ol.2 



l + dS - 92 1 + /5 1 - 5 ' 



(2) 



(3) 



(4) 



Here, Ko is the Kolmogorov constant (assumed to be 1.8), 
S — 162 A 2 (V — V a( j), and the coefficients are given by 
a = 10.8654, b = 4.89073xl0~ 3 , k = 0.149888, I = 0.189238, 
q = 1.85011, c = 1.08071 x 10~ 2 , d = 3.01208 x 10~ 3 , e = 
3.34441 x 10~ 4 , and / = 1.25 x 10~ 4 . A is given by 



c p p 2 KZ 2 ( g8 



12acT 3 ^ 2H D ^ 



1/2 



(5) 



A important characteristic of the CGM model is the 
absence of free parameters. In particular, the mixing length 
is taken as I = z where z is the geometrical distance from 
the top of the convection zone. Stothers & Chin (1997) have 
recently shown that I — z performs very well in different kind 
of stars. In order to avoid numerical difficulties, we allowed 
for a very small overshooting at the moment of evaluating 
I, that is, we considered I = z + /3Hl° p , where Jf* op is the 
pressure scaleheight at the top of the outer convection zone. 
We elected /3 < 0.1, which does not affect the temperature 
profile of the convection zone of our evolving models. 

For the sake of comparison, we also included in our cal- 
culations the ML1, ML2 and ML3 parametrizations of MLT. 
These versions, which are associated with different convec- 
tive efficiencies, has been amply employed in white dwarf 
investigations (see Tassoul at al. 1990 for details). 
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3 EVOLUTIONARY RESULTS 

Here, we present the main results of our calculations. We 
computed the evolution of DA white dwarf models with 
masses ranging from M = 0.5Mq to M — I.OMq at in- 
tervals of O.IMq and with a metallicity of Z — 0.001. We 
varied M H within the range 10" 13 < M H /M < 2 x 10" 4 
and Mho in the range 10 -6 < M He /M < 10" 2 . We employ 
the FST approach given by CGM, and also the ML1, ML2 
and ML3 versions of MLT. The models were evolved from 
the hot white dwarf stage down to log (L/Lq) = —5. Un- 
less stated otherwise, we shall refer in what follows to those 
models parametrized by M H = 10" 6 M and M Hc = 10~ 2 M. 

We begin by examining Fig. 1 in which the surface 
gravity of our DA models with logg(H) = —10,-6,-4 
and no hydrogen envelope is depicted in terms of the ef- 
fective temperature (q is the external mass fraction given 
by q = 1 - M r /Mf\ The first observation we can make 
from this figure is that thick hydrogen envelopes apprecia- 
bly modify the surface gravity of no - hydrogen models, es- 
pecially in the case of low - mass configurations. The effect 
of adding an outer hydrogen envelope of AIh = 10 -10 M 
is barely noticeable. Note also that finite - temperature ef- 
fects are relevant even in massive models (see also Koester 
& Schonberner 1986). At low effective temperatures, convec- 
tive mixing between hydrogen and helium layers (see later 
in this section) change the surface gravity of models with 
thin hydrogen envelopes (see the discussion in Benvenuto 
& Althaus 1998 in the context of low - mass helium white 
dwarfs for more details about this topic). 

It is worthwhile to mention that our more massive mod- 
els begin to develop a crystalline core at effective temper- 
atures close to the temperature range of the ZZ Ceti in- 
stability strip. More specifically, the onset of crystallization 
at the centre of 1.0, 0.9, 0.8 and 0.7 Mq DA models with 
log 5 = -6 takes place at T cff = 16200, 13100, 10870 and 
8900 K, respectively. For details concerning the process of 
growth of the crystal phase of our models, see Benvenuto & 
Althaus (1995). 

The role played by hydrogen burning in our evolving 
models is worthy of comment. Iben & Tutukov (1984) were 
the first in drawing the attention to the fact that hydrogen 
burning in cooling white dwarfs could be, at low luminosi- 
ties, a major energy source. We find in particular that for 
models with logg(H) > —4, hydrogen burning may con- 
tribute non - negligibly to the surface luminosity. For more 
details, we depict in Fig. 2 the relative contribution of hydro- 
gen burning as a function of effective temperature for various 
stellar masses having a hydrogen envelope of log g(H) = —4. 
Note that for intermediate stellar masses the hydrogen - 
burning contribution reaches a maximum near the ZZ Ceti 
instability strip. Another feature worthy of mention shown 
by Fig. 2 is that, even for more massive models, hydrogen 
burning never becomes the dominant source of stellar en- 
ergy. We also include in Fig. 2 the neutrino luminosity of 
the models. For the effective temperature range shown in the 
figure, neutrino luminosity strongly decreases with cooling 



T For the sake of clarity, in this and the following figures, we do 
not depict the results corresponding to 1.0 Mq models 
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Figure 1. Surface gravity versus the effective temperature for 
DA white dwarf models with (from top to bottom) M/Mq = 0.8, 
0.7, 0.6 and 0.5. For each stellar mass, the four curves correspond 
to models with (from top to bottom): no hydrogen envelope, and 
logg(H)= -10, -6 and -4. Note that very thick hydrogen envelopes 
appreciably modify the surface gravity of the no - hydrogen mod- 
els. 
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Figure 2. Hydrogen - burning (solid lines) and neutrino lumi- 
nosity (dotted lines), (in terms of the surface photon luminosity) 
as a function of the effective temperature for 0.5, 0.6, 0.7, 0.8 and 
0.9 Mq DA white dwarf models with log q(U) = -4. 



and it is negligible by the time models enter the instablitiy 
strip. 

The importance of nuclear burning is strongly sensitive 
to the exeat value of the hydrogen layer mass, as was recog- 
nised by Koester & Schonberner (1986) (see also D'Antona 
& Mazzitelli 1979). To show this, we have computed addi- 
tional sequences by considering the cases with logq(H) = 
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Figure 3. Hydrogen - burning luminosity (in terms of the surface 
photon luminosity) as a function of the effective temperature for 
0.6 and 0.8 Mq DA white dwarf models. For each stellar mass 
and from top to bottom, models with logg(H) = -3.699, -3.824 
and -4 are depicted. 



Figure 4. The location of the top and the base of the outer 
convection zone expressed in terms of the mass fraction q versus 
the effective temperature for a 0.6 Mq DA white dwarf model 
according to different theories of convection. For the small jump 
at T eS Ri 9000 K see text. 



-3.824; -3.699. The results for 0.6 and 0.8 M models are 
illustrated in Fig. 3. It is clear that a small difference in 
the hydrogen layer mass is responsible for the different role 
of hydrogen burning. In particular, for the 0.6 Mq model 
with Mh = 9 x 10~ 5 Mq, we find that the relative con- 
tribution of hydrogen burning at low luminosities remains 
always below 9 per cent, while for the same model but with 
Mh = 1.2 x 10~ 4 Mq, the hydrogen - burning contribution 
rises up to « 18 per cent 

The behaviour of the evolving outer convection zone 
is depicted in Figs. 4-7 for DA white dwarf models with 
M = 0.5, 0.6, 0.8 and 1.0 Mq and for the different the- 
ories of convection. In each figure, we plot the location of 
the top and the base of the outer convection zone in terms 
of the mass fraction q as a function of effective tempera- 
ture. The first observation we can make from these figures 
is that at both hot and cool extremes the convection zone 
profile is independent of the treatment of convection. In fact, 
at high effective temperatures only a negligible fraction of 
energy is carried by convection and, consequently, the tem- 
perature stratification remains essentially radiative, while 
at low temperatures convection is very efficient and most of 
the convection zone assumes an adiabatic stratification. At 
intermediate effective temperatures, however, where the ZZ 
Ceti stars are observed, convection treatment is decisive in 
fixing the structure of the outer zone of the models. Note 
that, in contrast to the situation encountered for DB ob- 
jects (Benvenuto & Althaus 1997), the CGM predictions are 
roughly intermediate to those of the ML1 and ML2 convec- 
tion. Another feature shown by these figures and worthy of 
comment is that the thickness of the convection zone in the 
CGM model begins to increase at a given effective temper- 
ature much more steeply than in any of the MLT versions. 
A similar trend is also found in DB white dwarf models, 
though it is worth remarking that the global behaviour of 




16 14 12 10 8 6 4 



T eff [x 10 K] 

Figure 5. Same as Fig. 4, but for a 0.5 Mq DA white dwarf 
model. 

the evolving convection zone is markedly different in both 
types of stars. In agreement with previous results (D'Antona 
& Mazzitelli 1979), the final extent reached by the base of 
the convection zone is smaller in the more massive mod- 
els. We should mention that the small jump at the top of 
the convection zone at T e g ~ is brought about, not 

surprisingly, by the change from OPAL to Alexander & Fer- 
guson's (1994) molecular opacities. This jump affects the z 
values, thus giving rise to some irregularities in the CGM 
convection profile, which does not alter the conclusions of 
this work. 
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Figure 6. Same as Fig. 4, but for a 0.8 Mq DA white dwarf 
model. 
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Figure 7. Same as Fig. 4, but for a 1.0 Mq DA white dwarf 
model. 



To clarify better the role played by convective efficiency 
in determining the thermal structure of the outer layers, wc 
show first of all in Fig. 8 the dimensionless convective tem- 
perature gradient in terms of the mass fraction q for the 
outer layers of a 0.6 Mq DA model with T eS = 10890^ 
according to the CGM model and the ML1 and ML2 ver- 
sions of MLT. Note the presence of a strong peak in V at 
logg ~ —16.1 in the CGM case. This arises simply because 
z < H p in the outermost layers, thus giving rise to very 
inefficient convection. In somewhat deeper layers, however, 
the CGM model turns out to be more efficient than ML1. 
This is a result of the fact that the larger number of ed- 
dies characterizing the CGM model, compared with MLT, 
begins to play a significant role. For \ogq > —15.6, both 
the CGM and ML2 models provide an adiabatic stratifica- 



Figure 8. The convective temperature gradient as a function 
of the mass fraction q for the 0.6 Mq DA white dwarf model at 
T cff = 10890-fiT according to the CGM, and to the ML1 and ML2 
versions of the MLT. The fact that I < H p in the outermost layers 
leads to a strong peak in V in the case of the CGM model. 




16.4 16.0 15.6 15.2 14.8 14.4 14.0 
-Log(q) 

Figure 9. The outer layer temperature versus the mass fraction 
q. The DA model corresponds to that of Fig. 8. Note that the 
temperature profile given by the CGM cannot be reproduced with 
an MLT with a single value of a. 



tion and hence the behaviour of V in the outermost layers 
will be critical for determining the temperaute profile in the 
envelope and thus the extent of the convection zone of the 
model. The resulting temperature profile of the model anal- 
ysed in Fig. 8 is depicted in Fig. 9. It is clear that the CGM 
model predictions are very different from those of the ML1 
and ML2 models. Indeed, the temperature profile cannot be 
reproduced with an MLT model with a single value of a. 

A final remark about the convection profiles is that at 
low effective temperatures they ultimately reach a maximum 
depth, which is determined by the location of the boundary 
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Table 1. Mixing temperatures and final hydrogen abundance for different models and convection theories 



Sequence 
0.50 (-12) 
0.50 (-10) 
0.50 (-8) 
0.50 (-6) 
0.60 (-13) 
0.60 (-12) 
0.60 (-11) 
0.60 (-10) 
0.60 (-8) 
0.60 (-6) 
0.70 (-12) 
0.70 (-10) 
0.70 (-8) 
0.80 (-12) 
0.80 (-10) 
0.80 (-8) 
0.90 (-12) 
0.90 (-10) 
0.90 (-8) 
1.0 (-12) 
1.0 (-10) 
1.0 (-8) 



ML1 
10,000 (<lc-6) 

9130 (le-5) 

7090 (2c-3) 
5600 («0.55) 
10,520 (<lc-6) 
10,130 (<lc-6) 

9680 (3e-6) 

8950 (3e-5) 
6710 (7.3c-3) 

5330 («0.8) 
10,160 (lc-6) 

8980 (8e-5) 
6410 (0.027) 
10,250 (2c-6) 
8740 (2.2e-4) 

6320 (0.1) 
10,350 (4c-6) 
8560 (6.5e-4) 

6170 («0.4) 
10,290 (le-5) 

8240 (2e-3) 

5950 (^0.8) 



ML2 

11,240 (<le-6) 
10,140 (9e-6) 
7450 (1.7e-3) 
5605 («0.55) 
11,820 (<lc-6) 
11,380 (<le-6) 
10,840 (2e-6) 
10,090 (2e-5) 
6900 (6.4e-3) 
5360 («0.8) 
11,390 (lc-6) 
9890 (6c-5) 
6530 (0.025) 
11,500 (lc-6) 
9740 (1.7c-4) 
6330 (0.1) 
11,470 (3e-6) 
9320 (5e-4) 
6200 (ra0.4) 
11,560 (7c-6) 
8810 (1.6e-3) 
5970 (^0.8) 



ML3 

11,870 (<lc-6) 
10,700 (9c-6) 
7470 (1.7c-3) 
5610 («0.55) 
12,505 (<lc-6) 
12,000 (<le-6) 
11,450 (2e-6) 
10,650 (2c-5) 
6850 (6.8c-3) 
5380 («0.8) 
12,100 (lc-6) 
10,330 (6c-5) 
6530 (0.026) 
12,020 (lc-6) 
10,010 (1.6c-4) 
6350 (0.1) 
12,050 (3c-6) 
9490 (4.6e-4) 
6220 («0.4) 
11,980 (7c-6) 
8950 (1.6c-3) 
5990 (f«0.8) 



CGM 
10,400 (<lc-6) 

9500 (le-5) 
7220 (1.7c-3) 
5610 (ra0.55) 
10,910 (<lc-6) 
10,570 (<lc-6) 
10,020 (3e-6) 

9450 (2e-5) 
6780 (6.6e-3) 

5360 (ra0.8) 
10,640 (lc-6) 

9320 (7c-5) 
6500 (0.025) 
10,760 (lc-6) 
9180 (1.8e-4) 

6350 (0.1) 
10,830 (3e-6) 

9120 (5e-4) 

6230 (psO.4) 
10,800 (8e-6) 
8610 (1.6e-3) 

6000 («0.8) 



Each sequence is denoted by the stellar mass (in solar mass units) and the value of log qjj . The numbers in 
parentheses next to each temperature give the hydrogen surface content after mixing. 



of the degenerate core. For the 0.60 Mq DA white dwarf 
model, in particular, such a depth occurs around log q ft* —6. 
(see Tassoul et al. 1990 for a similar result). Accordingly, 
mixing episodes between hydrogen - rich and the underlying 
helium - rich zones are possible only for values of log g(H) 
lower than this limit. It is apparent from Figs. 4-7 that the 
mixing effective temperature for models with very thin hy- 
drogen envelopes will be strongly dependent on the assumed 
convective efficiency and not so much on the stellar mass. 

To put this assertion on a more quantitative basis, we 
list in Table 1 the effective temperatures at which the hydro- 
gen surface abundance begins to decrease appreciably as a 
result of convective mixing, together with the final hydrogen 
abundances (for models with thick hydrogen layers, mixing 
process takes place gradually in a finite range of effective 
temperature; in that case we only give approximate values 
of the hydrogen abundance). Clearly, for models with very 
thin hydrogen envelopes, convective mixing drastically mod- 
ifies the surface composition from a hydrogen - dominated 
to a helium - dominated one. This behaviour can be under- 
stood directly by examinig Fig. 10, which shows the evolving 
convection zones of the 0.6Mq model with log q(R) = —12 
according to the CGM convection. At T e g « 10600AT , the 
merging of the hydrogen convection zone with the helium 
convection zone (located just below the hydrogen/helium 
transition region) causes the base of the convection zone to 
reach very deep helium layers quickly, giving rise to an al- 
most complete dilution of the hydrogen content. From then 
on, the subsequent evolution corresponds to that of a DB 
model (see Fig. 11 of Benvenuto & Althaus 1997). 

In order to estimate approximate effective temperatures 
for the theoretical blue edge of DA instability strip we em- 
ploy thermal time - scale arguments 7th for our evolving 
models. Several past studies (Cox 1980, Winget et al. 1982, 
Tassoul et al. 1990 and references cited therein) have shown 
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Figure 10. The location of the top and the base of the outer 
convection zones expressed in terms of the mass fraction q versus 
the effective temperature. The results correspond to a 0.6 Mq 
DA white dwarf model with logq(H) = —12 according to the 
CGM model. Note the presence of two convection zones at high 
T e g which ultimately merge at T ~ 10500-K", thus giving rise 
to a drastic change of the surface chemical composition from a 
hydrogen - dominated to a helium - dominated one. 



that the behaviour of pulsational instabilities in white dwarf 
stars can be understood in terms of r th of the driving regions 
defined by 



Tth 



JO 



CV T 



M* dq 



(6) 
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Table 2. Theoretical blue-edge effective temperatures versus stel- 
lar mass and convection theory for our DA white dwarf 

Theory of convection Mass (M/Mq ) T e ff (K) 
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In Eq. 6, Cy is the specific heat at constant volume, L 
and M are, respectively, the luminosity and the mass of the 
model, and the subscript "be" corresponds 

to the location of the bottom of the superficial convec- 
tion zone. In particular, the blue edge of the instability strip 
(the onset of pulsations) corresponds approximately to the 
effective temperature at which r t h becomes comparable to 
100 s, which for a white dwarf is of the order of the shortest 
observable g - mode periods. In this connection, we list in 
Table 2 (see also Fig. 11) the effective temperature of our 
theoretical blue edges of the ZZ Ceti instability strip ac- 
cording to the different stellar masses and convection treat- 
ments we have considered. Note, as found by other investi- 
gators (e.g. Tassoul et al. 1990 and Bradley & Winget 1994) 
the dependence of the blue edge temperature on the stel- 
lar mass and, more importantly, on the assumed convective 
efficiency. In particular, our MLT blue edges are consistent 
with those of Bradley & Winget (1994) obtained on the basis 
of detailed pulsation calculations. Clearly the CGM theory 
predictions are intermediate to the ML1 and ML2 results. 
In Fig. 11, we have also included the theoretical result ob- 
tained by Gautschy et al. (1996) for a 0.6Mq ZZ Ceti model 
and the observational data of Bergeron et al. (1995) as well 
(we picked out those ZZ Ceti stars analysed by Bergeron et 
al. which, for a given stellar mass, are characterized by the 
highest temperature). 

With regard to the inclusion of the observational data 
in this figure, some comments seem to be appropriate at this 
point. As mentioned earlier, the effective temperature of the 
stars that define the observational blue edge is computed, 
as usual, by employing the emergent spectrum of a model 
atmosphere, which in the case of DA white dwarfs is strongly 
dependent on the prescription adopted for the treatment of 
convection (see e.g. Bergeron et al. 1995). Thus, in order to 
perform a self - consistent comparison with observations, we 



T 
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14 13 12 11 10 9 

T eff [x 10 3 K] 

Figure 11. The dependence of theoretical blue edge temperature 
on the stellar mass for CGM, and for the ML1, ML2 and ML3 
versions of the MLT. The filled triangles represent the observa- 
tions of Bergeron et al. (1995) for (from top to bottom) G226-29, 
G185-32, R548 and G238-53. The filled square corresponds to the 
blue edge temperature obtained by Gautschy et al. (1996) from 
hydrodynamical simulations of convection. Note that both the 
CGM model and the hydrodynamical simulation tend to point 
towards a cooler blue edge than the observed one. 

should employ a model atmosphere computed considering 
the CGM theory. In view of the fact that such models are 
still not available, any comparison of our theoretical results 
with observation should be taken with caution. 

In spite of the above warnings, it is remarkable that 
both the CGM theory and Gautschy et al.'s (1996) predic- 
tions tend to point towards a cooler blue edge than obser- 
vations. In particular, on the basis of very detailed hydro- 
dynamical simulations of convection in a 0.6Mq ZZ Ceti 
model, Gautschy et al. derived an effective temperature be- 
tween 11400 K and 11800 K for the blue edge. As far as the 
blue edge dependence on the stellar mass is concerned, the 
MLT and CGM theory predictions are rather similar, span- 
ning « 1QQQK over the range of stellar masses we considered. 
Note that dependence of the blue edge on the stellar mass 
predicted by both theories of convection is rather similar to 
that shown by observations. 



4 SUMMARY 

In this paper we compute the structure and evolution of car- 
bon - oxygen white dwarf models with hydrogen envelopes 
(DA type). We consider stellar masses ranging from M=0.5 
to M — 1.0 Mq at intervals of M = O.IMq, and we treat 
the mass of the hydrogen and helium envelopes as free pa- 
rameters within the range 10" 13 < M H /M < 10~ 4 and 
10~ 6 < Mho/M < 10~ 2 , respectively. The models were 
evolved from the intermediate effective temperature stage 
down to log (L/Lq) — —5. 

The calculations were made with a white dwarf evo- 
lutionary code including updated radiative and conductive 
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opacities and neutrino emission rates, and very detailed 
equations of state for hydrogen and helium plasmas. We 
also include the effects of crystallization, convective mix- 
ing and hydrogen burning (pp chain and CNO bi - cycle). 
The most important feature of our study, however, is that 
we treat the energy transport by convection within the for- 
malism of the so - called full - spectrum turbulence theory, 
which constitutes an improvement over most previous stud- 
ies of DA white dwarf evolution. In particular, we employ 
a new model based on such theory developed by Canuto, 
Goldman & Mazzitelli (1996) (CGM). This new model in- 
cludes the full spectrum of eddies, has no free parameter 
and computes the rate of energy input self - consistently. In 
the white dwarf domain, the CGM models has been recently 
shown (Althaus & Benvenuto 1997b) to provide a good fit to 
new observational data of pulsating DB objects. Finally, for 
the sake of comparison, we also consider the most common 
parametrizations of the mixing length theory (MLT) usually 
employed in this kind of studies. 

In agreement with previous studies, we find that very 
thick hydrogen layers substantially modify the surface grav- 
ity of the models and that the importance of nuclear burning 
is strongly sensitive to the exeat value of the hydrogen layer 
mass. In particular, we find that for the 0.6 Mq model with 
Mh = 9 x 10~ 5 Mq the relative contribution of hydrogen 
burning at low luminosities remains always below 9 per cent, 
while for the same model but with Mh = 1.2 x 10 -4 Mq, the 
hydrogen burning contribution rises up to ~ 18 per cent. 

One of our main interest in this work has been to study 
the evolution of ZZ Ceti models with the aim of comparing 
the CGM and MLT predictions. In this connection, we find 
that the temperature profile given by the CGM model is 
markedly different from that of the ML1 and ML2 models, 
and it cannot be reproduced by any MLT model with a single 
value of a. The evolving outer convection zone also behaves 
differently in both theories. In particular, the thickness of 
the convection zone in the CGM model begins to increase 
at a given effective temperature much more steeply than in 
any of the MLT versions and remains intermediate to those 
with ML1 and ML2 convection. 

We have also computed approximate effective tempera- 
tures for the theoretical blue edge of DA instability strip by 
using thermal time - scale arguments for our evolving DA 
models. In this context, we find that the effective tempera- 
ture of the theoretical blue edges of the ZZ Ceti instability 
strip depends strongly on the stellar mass and more impor- 
tantly on the assumed convective efficiency, in agreement 
with previuos studies. In particular, our MLT blue edges 
are consistent with those of Bradley & Winget (1994) ob- 
tained on the basis of detailed pulsation calculations. We 
find that the CGM theory predicts blue edges cooler (by 
~ lOOOA') than the observed ones. Specifically, we find an ef- 
fective temperature of llOOOif for the blue edge of a 0.6 Mq 
DA model. It is remarkable that recent non - adiabatic pul- 
sation calculations based on numerical simulations of con- 
vection tend also to indicate a somewhat cooler blue edge. 
However, we remaind the reader that the physical ingredi- 
ents we employed in the present paper are not consistent 
with those of the stellar atmosphere calculations employed 
by Bergeron et al. (1995), particularly in the treatment of 
convection. Although a difference of 1000 K between the 
effective temperature of the observed and theoretical blue 



edges for the ZZ Ceti instability strip seems to be large, 
we should wait for model atmospheres also computed in the 
frame of the full spectrum turbulence theory. Only after such 
models become available we shall be able to gauge the actual 
importance of such a discrepancy. 

Detailed tabulations of the evolution of our DA mod- 
els, which are not reproduce d here, are available at 
http:/ /www.fcaglp.unlp.edu.ar/ ~althaus/ 
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